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ABSTRACT 

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from 
the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI 
DR1). The clouds were identified by a custom machine- vision algorithm utilizing 
Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds 
have velocities typically between |Vlsr| = 20 — 400 km s -1 , linewidths of 2.5 — 35 
km s _1 , and column densities ranging from 1 — 35 x 10 18 cm -2 . The distances to 
the clouds in this catalog may cover several orders of magnitude, so the masses 
may range from less than a Solar mass for clouds within the Galactic disc, to 
greater than 10 4 M for HVCs at the tip of the Magellanic Stream. To search for 
trends, we separate the catalog into five populations based on position, velocity, 
and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity 
clouds (LVCs); warm, low positive- velocity clouds in the third Galactic Quadrant; 
and the remaining warm LVCs. The observed HVCs are found to be associated 
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with previously- identified HVC complexes. We do not observe a large population 
of isolated clouds at high velocities as some models predict. We see evidence for 
distinct histories at low velocities in detecting populations of clouds corotating 
with the Galactic disc and a set of clouds that is not corotating. 

Subject headings: Galaxy: halo — intergalactic medium — galaxies: formation — 
cooling flows — Galaxy: ISM 



1. Introduction 



The recent discovery of compact clouds over the entire velocity range of our Galaxy 



tion and the Galactic 


ISM (e.g. 


Lockman 


2002 


>■ 


Ford et al. 


2008 




Rvan- Weber et al. 


2008; 


Ben Bekhti et a 


J 2009 


Stanimirovic et al. 


2008; 


Heitsch & Putman 


20091: 


Begum et al. 


2010; 


Giovanelli et al. 


2010; 


ilsu et al. 


2011 


). Despite the similarity of the individual clouds when 



observed with the hyperfine neutral hydrogen 21 cm line, different classes of clouds have 
been identified and many possible origins have been proposed. The detection of stellar out- 
flows in H I points to the possibi lity that some small clouds could be produced by evolved 



stars (e.g. iMatthews et al.ll201ll ). At the disc-halo interface, a population of clouds has 



been observe d rotating with the Gal axy that may have been pushed off the disc by stel- 



lar feedback (jFord et al 



2008 



20101 ). In the Galactic halo, compact clouds may repre- 
sent the initial cooling seeds from the multiphase halo medium that will eventually fuel 
the disk (IMaller fc Bullock! 12004 ; iJoung et al.l 120111 ). or the remnants of larger complexes 
( Heitsch fc Putman 12009 ). The continual production of these clouds could provide a sig- 
nificant source of the fuel for the Milky Way's ongoi ng star formation. The abundance 



of newly discovered sm all Local Group dwarf galaxies ( IBelokurov et al.l l2Q06t IZucker et al. 



2006 



Irwin et al.ll2007l . and others), one of which is particularly gas- rich (IRyan- Weber et al. 
20081 ) , raises the possibility that some of these clouds represent previously undiscovered dwarf 
galaxi es. It is even possible that some compact clouds are outside the halo and physically 
large ( IGiovanelli et al.l 120101 ). 



Unraveling the nature of the compact clouds requires a large, well-defined sample. His- 
torically, finding such a sample has been difficult due to the limited area, sensitivity, and 
resolution covered by HI surveys, a s well as the difficulty of recoverin g compact objects con- 
fused by diffuse Galactic emission (jStil et al.ll2006al : iFord et al.ll2008l ). The latter problem is 
evident in the lack of detections of halo clo uds (commonly referred to as high- velocity clouds ; 
HVCs) between V LSR = -90 to 90 km s' ^Wakker fc van Woerden1ll997l : Ipeek et al.l l2009h . 
and the recent, optical discovery of gas-rich dwarf satellites at velocities where Galactic 
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emission is prevalent ( llrwin et al.l 120071 ). A large, sensitive survey with sufficient kinematic 
and spatial resolution to recognize clouds that are near Galactic emission is needed. 



The GALFA-HI (Galactic Arecibo L-Band Feed Array HI) Survey flPeek et al.l 120111 ) 
provides the data needed to detect a large sample of compact clouds across the range of 
velocities of the Galaxy and Local Group. At completion it will cover 13,000 deg 2 (58% 
of which is available in Data Release iQand cataloged here; DR1) over -1°< 5 < 38° and 
the entire RA range between Vlsr= ±650 km s _1 . The GALFA-HI survey has a spatial 
resolution of 4' and a channel spacing of 0.184 km s _1 . The sensitivity at 6a to compact, 
warm (FWHM = 15 km s _1 ) HI clouds corresponds to masses as low as 7 x 10~ 2 M 
at 1 kpc and 7 x 10 4 M at 1 Mpc (see §2]). As the sensitivity of a survey to compact, 
spectrally resolved sources scales inversely with both the antenna beam solid angle and the 
noise for a given channel width, GALFA-HI is uniquely positioned among other large-area 
Galactic HI surveys. With the sensitivity of a single-dish radio telescope and resolution 
approaching that of a compact array, GALFA-HI is 2. 5 to 80 times more sensitive to com 



pact clouds than other large- area surveys (2.5: IGPS: iTavlor et al. 



McClure- Griffiths et al. 



2010 



2005 



2003 



15: GASS / EBHIS: lMcClure-Griffiths et al. 



Stil et al. 



2009 



2006b 



Winkel et al. 



80: LAB: Kalberla et al.l 120051 ). 



through an automated algorithm (IHsu et al.ll201l[ ) 



nitial investigations of GALFA-HI data id enti 



fied many compact clouds, both v isually (IBegum et al.ll2010l : IStanimirovic et al.l 120061 ) and 



In this paper we present a catalog of compact clouds (4'-20') from the GALFA-HI 
Survey DR1 that has been created using a new technique designed to recover faint clouds 
and identify structures near diffuse Galactic emission. The methods we use are detailed in 
§|3l with the completeness of the catalog assessed in §3.31 and the catalog presented in §3.41 
The properties of the clouds cataloged are described in §@] where the clouds are separated 
into five populations. A discussion of the implications of our results can be found in §5J 



2. Observations and Data Reduction 



The catalog is generated using the GALFA-HI Survey DR1. GALFA-HI is a survey of 
the 1420 MHz hyperfine transition of neutral hydrogen in the Galaxy using the Arecibo 305m 
telescope and the ALFA seven-beam feed array. Th e survey is completed c ommensally with 
other Arecibo extragalactic and Galactic surveys ( jGiovanelli et al.l l2005t iGuram fc Taylor 
20091 ). GALFA-HI data provide a channel spacing of 0.184 km s _1 and cover a velocity range 
of ±650 km s _1 in the Local Standard of Rest (LSR) with a spatial resolution of 4'. The DR1 



lr The GALFA-HI DR1 is publicly available at 



http://purcell.ssl.berekley.edu 
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data cover 7520 square degrees of sky in an area between 5 = 38° and 5 = — 1° (see bottom 
panel of Figure HJ and Figure |2]), with a range of sensitivity from 120mK to 50mK in 0.74 
km s _1 channels. The data were corrected for the effects of bandpass, gain, baseline ripple, 
and the first sidelobes of the ALFA beams with a special-purpose data reduction pipeline. 
The data are not stray radiation corrected. As the Arecibo stray sidelobes are large and low 
gain, the effect of stray radiation is to a dd a spurious low -amplitude HI signal that varies 
slowly as a function of position on the sky (IPeek et al.ll20!ll ). Were A recibo to have angularly 
compact stray sidelobes, broadband surveys such as GALFACTS (IGuram fc Taylorl 120091 ) 
would detect erroneous replications of tremendously bright continuum sources, an effect 
that has not been obeserved. The spatially broad contamination from stray radiation should 
not significantly affect our catalog of compact clouds. The details of GALFA-H I observing 



and d ata reduction, along with the specifics of the DR1 data set, can be found in lPeek et al. 



3. Cataloging Method 

Identifying clouds by eye, as has been done in the past, is impractical for the amount 
of data the GALFA-HI survey is producing. This, coupled with a desire for completeness 
over a broad range of velocities and sensitivities, motivated our development of an algorithm 
to identify and characterize compact isolated clouds in the GALFA-HI data cubes. The 
algorithm consists of four distinct components: Galactic Subtraction, Convolution, Region 
of Interest (ROI) Detection, and Merge and Candidate Selection; each discussed in detail in 
§3.11 After the initial detection, each cloud candidate was inspected to remove systematic 
noise spikes, RFI, and connected structures. This parsing process is discussed in §3.21 In 
§3.31 we show that our sensitivity is dominated by systematic errors and that our sensitivity 
is effectively constant over our search space. The properties we measure for each cloud are 
explained in §3.41 



3.1. Cloud Detection Algorithm 

1. Galactic Subtraction - To remove smooth Galactic emission, we applied a high-pass 
spatial filter by smoothing each channel with a one square degree median box and subtracting 
the resulting smoothed cube from the data. The intent of this pre-processing step is to 
remove emission on larger scales than our targets (20') without removing signal from the 
compact clouds. We chose a box size of one degree because an isolated 20'cloud on a smooth 
background is unaffected by the filtering. We applied the filter over the full velocity range 
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to avoid introducing discontinuities in the data. Filtering allows us to probe for compact 
clouds down to velocities closer to Vlsr= km s _1 than would otherwise be possible. We 
do not make a lower cut in velocity. Instead, we found the amount subtracted in this step 
to be a good metric for cloud belie vability. We discuss this further in §3.21 

2. Convolution - To increase our sensitivity, we convolved the data with a three dimensional 
kernel in a simplified wavelet analysis. Using this method enhances regions of isolated 
emission according to their similarity to the 3D kernels used. We chose to use four "Difference 
of Gaussians" (DOG) kernels, one of which is shown in Figure [3J The DOG kernel acts like 
a bandpass filter and is computationally fa ster than other ker nels since it can be separated 



into ID and 2D kernels used in series (See ISonka et al.l 120071 . and references therein). The 



kernel is the sum of two three-dimensional Gaussians: a positive Gaussian and a negative 
Gaussian that is wider in all dimensions by a factor of 1.2. The Gaussians are scaled such 
that the integral of the kernel is zero. We define the size of the kernel as the FWHM of the 
positive Gaussian. 

We tuned the sizes and widths of the four kernels by measuring the peak response 
to a set of model clouds of fixed signal to noise. The model clouds covered the range of 
possible observed cloud sizes and velocity widths given the data and chosen upper limits. 
The response of the optimum kernels is plotted in Figure HJ The four kernels are circular 
with FWHM and velocity widths of 7' and 5 km s _1 , 7' and 15 km s _1 , 18' and 5 km s _1 , 
and 18' and 15 km s" 1 . The theoretical response is above 80% for our entire search space 
except for clouds that are both large (20') and wide (20 km s _1 ) where the response drops to 
70%. While increasing the number of kernels would provide a smoother sensitivity function, 
the computational requirements are unreasonable. §3.31 discusses sensitivity further. 

This step results in four convolved data cubes that will be passed to the region detecting 
routines. Before this, a careful treatment of noise is necessary to accurately identify levels 
of significance and to make it possible to compare the results for the four kernels. The noise 
characteristics are complicated due to different regions of D R1 having varying amounts of 



data taken with different observing modes (see Figure [Tl and iPeek et al.l (l201lh for details). 
Also, the noise in each voxel is not statistically independent. While the GALFA-HI data 
cubes are gridded to 1' pixels, Arecibo's resolution is 4'. We also observe correlated noise 
on the scale of ALFA's beam pattern of ~25'. In order to calculate the uncertainty in the 
convolved cubes we take the following steps. (1) We begin by taking the original data and 
finding the standard deviation along each line of sight over a velocity range that does not 
include significant emission. This allows us to construct a noise map for the original data. 
(2) Assuming the the value of the noise is the same over the full velocity range for each pixel, 
the noise detected by each kernel is calculated by multiplying the noise map by each kernel 
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at each position, and then summing the scaled expected noise values in quadrature. This 
creates a convolved noise map for each kernel. (3) Dividing each channel in the convolved 
cubes by the corresponding convolved noise map converts the cube to units of signal to noise, 
for uncorrelated noise. (4) We normalize the data to correct for the effect of correlated noise 
in the unconvolved datacubes. We fit a Gaussian to the distribution of values in a region 
with no detectable emission. For uncorrelated noise, the standard deviation should be one. 
Dividing each convolved cube with the corresponding fit's standard deviation converts the 
units to signal to noise (a). 

3. Region of Interest (ROI) Detection - We now have four convolved data cubes in 
units of a that can be used to detect regions of interest (ROI). We use a simple 'watertable' 
technique in which the highest peak in a cube is identified and the adjacent voxels above a 
watertable of 4<r are added to create the ROI. The ROI is then numbered, saved, and deleted 
from the cube. This is then repeated with the new highest peak in the cube selected and 
the ROI grown to the watertable level again. The ROI detection process is completed when 
there are no peaks above a minimum of 6a. 

4. Merge and Candidate Selection- Having defined ROIs in all cubes corresponding to 
the four kernels, we need to merge the information to select clouds. For this catalog we have 
chosen to include clouds with a spatial FWHM less than 20'. Our first step is to reject those 
clouds larger than this selected by our kernels. For an ideal, Gaussian-shaped cloud with a 
spatial FWHM of 20' and a peak brightness of IK in the unconvolved data, the 7' kernels 
would identify a ROI 34' across and the 18' kernels would identify a ROI 36' across. Therefore 
ROIs larger than these values are rejected as extended, and any ROI's found within these 
larger structures in position- velocity space are flagged as "embedded" and also removed. We 
also removed ROIs that extend into velocities where the bandpass response begins to decline 
(|V LS r| > 650 km s" 1 ). 

We use four kernels of different sizes to achieve smooth sensitivity over the parameter 
space being searched. And while the weakest clouds are only detected by the most similar 
kernel, brighter clouds can be detected by multiple kernels. Therefore, the next part of 
the Merge and Candidate Selection step is to take the clouds found by multiple kernels 
and select the ROI for the cloud that has the highest peak value. A higher response from a 
specific kernel indicates a better fit to the cloud. The final list of ~10,000 ROIs were visually 
inspected and classified as described in the next section. 
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3.2. Cloud Classification 

The cloud detection algorithm is very effective at detecting significant isolated struc- 
tures, but the ISM includes extremely complicated features. At low velocities in particular, 
the Galactic emission is a web of needles, c ompact clouds at the en d of filaments (peninsu- 



las), bubbles, arcs, and many other shapes (IKalberla fc Kerpll2009l . and references therein). 



Some of these structures were included in our candidate list, as well as some RFI and cal- 
ibration artifacts. Removing connected structures and artifacts requires visual inspection 
of the candidates. When visual inspection began it was soon realized that those ROIs in 
regions where more than IK was removed in the Galactic Subtraction step were generally 
part of these more extended structures. We therefore rejected those ROIs, resulting in the 
gap in velocity coverage seen in Figures [fl El and|8] (generally |Vlsr| < 20 km s _1 ). 

The remaining ROIs were examined visually to determine if they were part of a more 
extended Galactic structure, RFI or a calibration artifact, and should be removed. To 
be included in the final catalog, we required that a cloud was visually inspected by at 
least three of the authors, with a majority in agreement. Since the clouds are detected 
by integrating three dimensional volumes, visual inspection is a non-trivial process. With 
thousands of clouds to inspect and ten viewers, bookkeeping was also non-trivial. To solve 
both of these problems, we built a procedure which displays integrated position-position 
images and average spectra similar to those in Figure If this information was not enough 
to classify the ROI, the user could scan through the channel images and spectra of the ROI 
in the original or Galactic-subtracted cubes. The viewer then classifies the ROI as a real 
discrete cloud, a cloud that is part of a more extended structure, or a glitch in the data. 
Only ROIs in the first category are included in the final cloud catalog. A ROI was considered 
part of a more extended structure and excluded from the catalog if there was emission > 4cr 
connecting the ROI to other areas of emission. This is 4er in the unconvolved cube and 
represents structures that are not well enough defined by the kernels to appear significant in 
the convolved cubes. If there are other areas of emission nearby, either in the form of other 
compact clouds or a larger complex, the cloud was included as long as it was not connected 
at the 4<t level. Clouds with average spectra so weak that the fitting described in §3.41 failed 
were also removed from the catalog. 

After pruning the catalog, the remaining clouds were compared to the NED database to 
remove any previously known galaxies. We searched for known objects within 50 km s" 1 in 
velocity and 10' in position for each of the catalog candidates. Those candidates in agreement 
with known galaxies were removed from the catalog. Of the 2025 candidate clouds after visual 
evaluation, 61 were found to be associated with known galaxies and were removed leaving a 
final catalog of 1964 clouds. 
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3.3. Empirical sensitivity 

To fully understand the cloud sample we must understand the limitations and sensitivity 
of the cataloging process. Since the background contaminants can deviate dramatically 
from Gaussian noise and our cataloging method is non-linear, it would be very difficult to 
accurately determine our sensitivity only using an analytic technique. As such, we test our 
sensitivity using a Monte Carlo method; we inject false clouds into the data set, run our 
entire cataloging process and find which of these clouds we can recover. The false clouds we 
inject are 3D Gaussians, with variable RA, Dec, velocity, peak brightness (in Kelvin; T p k), 
major axis (FWHM in arcminutes; a), aspect ratio (b/a), linewidth (FWHM in km s _1 ; AV) 
and position angle (in degrees; 0). We inject them into two regions, one with integration 
times (t) typically larger than 40 seconds per beam, and one with integration times typically 
smaller than 40 seconds per beam. We note here that many of the real clouds we detect in 
our catalog vary from simple Gaussian shapes, which will typically lower our sensitivity for 
a given set of parameters. 

We find that our detection fraction is entirely independent of position angle. In position 
and velocity, clouds are more difficult to detect when they are near the edge of the mapped 
area or in areas of strong emission where the Galactic Subtraction step removed more that IK 
(See §3.ip . Ignoring for the moment the effect of background HI on our sensitivity, we use the 
remaining parameters (T pk , a, b/a, AV, and t), to construct a metric, S (T pk , a, b/a, AV, t), 
to quantify the extent to which a cloud is detected. Since by construction we do not detect 
clouds whose value in the convolved cubes, <Jmax, is below 6a, we fit the functional form 
S (Tpk, a, b/a, AV, t) to match cr MAX , so as to best predict whether <tmax > 6a and thus be 
detected. We find that a simple exponential scaling, 

S {T pk , a, b/a, AV, t) = AT^(b/ af A\H e , (1) 

rather accurately captures the sensitivity function. We find A = 0.86a, a = 0.94, /3 = 
0.57, 7 = 0.18, 5 = 0.30, e = 0.50, using a least-squares fit, minimizing the difference between 
cmax and S. In Figure [6] we plot S (T pk , a, b/a, AV, t) / a max against each of the parameters, 
to demonstrate that there is very little functional dependence of S/ctmax to any of these 
parameters. We note that at values of a and AV below the size of our smallest kernel 
(5', 7 km s" 1 ) (Tmax decreases somewhat more quickly than the simple functional form in[IJ 
although the change is less than 50%. We find that we detect 90% of the clouds with S > 6a. 
The 10% we do not detect are typically coincident with other HI structures or glitches in the 
data, and thus not detected as compact. We found that the fraction of false clouds we expect 
to detect that go undetected is independent of background brightness, once we ignore clouds 
with background brightnesses above IK. This confirms that our background brightness cut 
of IK is relatively conservative. 
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Comparisons to previous catalogs pr ovide another metri c for our detecti o n algo rithm. 
We compare our catalog to the work of iBegum et al.l (120101 ) and iHsu et al.l (120111 ) . The 
clouds of the Begum et al. catalog were identified visually in the high-resolution GALFA-HI 
data cubes. They found 96 clouds with all but two at |Vlsr| < 90 km s -1 . We include 
31 of the Begum et al. clouds in our catalog. Of the remaining 65 clouds, 64 clouds were 
never selected as candidates due to being connected to extended emission. The one cloud 
that was selected as a candidate, but not included in our catalog, was removed because it 
was connected to the edge of an unobserved region. The median |Vlsr| of the undetected 
Begum et al. clouds is 31 km s _1 since many of these clouds are near Galactic emission. The 
positions, velocities, linewidths, and sizes of clouds in both catalogs are in agreement when 
the different methods used are considered. The measured velocities agree within 1 km s _1 and 
the typical size difference is less than 2'. Begum et al. use an average spectrum to measure 
brightness temperature and linewidth. They therefore measure lower peak brightnesses, and 
the typical difference in measured linewidth is less than 3 km s . 



We identify 26 of the Magellanic Stream clouds found by IHsu et al.l (120111 ). Of the 
remaining Hsu et al. Stream clouds, 19 are larger than 30' and are excluded, 15 were detected 
but were connected to other structure, and 15 were below our signal-to-noise threshold. The 
velocities and linewidths measured by Hsu et al. are tak en from an aver age spectrum over the 
region identified by the Duchamp detection algorithm (jWhitingll2008l ). Duchamp identifies 
regions over a given threshold (3<r for Hsu et al.) which are larger than the Aa regions we 
identify. This introduces differences in the measured velocities and line widths of typically 
5 km s _1 . The brightness temperatures are measured from the raw data and agree with 
our measurements to within 0.1 K. These results indicate that our algorithm is effective at 
identifying and characterizing clouds that meet our criteria: smaller than 20', brighter than 
6<r in the convolved data, and separated from complex, extended structure in position and 
velocity. 



3.4. Cloud Properties 

To characterize the 1964 identified clouds, we began by fitting a four parameter Gaus- 
sian to the average spectrum over the ROI. The four parameters were amplitude, velocity, 
linewidth, and a constant baseline. Using the results of the fit, we then refined the velocity 
range and fit the average spectrum over the new velocity range with the same single, four 
parameter Gaussian. The images in the catalog (See Figure [5] for an excerpt) were prepared 
by scaling each channel by the normalized average fit and then integrating over velocity. 
This produces cleaner images and deemphasizes nearby structures at velocities where the 
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clouds of interest are weakest. To calculate the maximum column density of the cloud, we 
fit the same type of single Gaussian to the location of the brightest peak in the integrated 
image. When the fit for the maximum column density fails we report the average column 
density lower limit. 

A sample of the catalog is presented in Table [1] and Figure [5j while the online version 
contains the complete catalog §The definitions of the cloud characteristics listed in Table [1] 
are summarized below. We note that the positions are obtained from the convolved cubes, 
while the other parameters are obtained from the original cubes after Galactic subtraction. 



• Cloud ID - Clouds are named by the right ascension (in degrees), declination, and 
Vlsr velocity of the peak in the convolved cube. 

• RA & DEC - The right ascension and declination of the peak in the convolved cube. 
Note that this is not necessarily the location of peak flux, or the spatial center of the 
region of interest (see §3. 1 1) . 

• I & b - The Galactic longitude and latitude corresponding to the right ascension and 
declination. 

• o"max - The maximum signal-to-noise of the cloud in the convolved cube. This is a 
measure of the significance of the cloud. 

• Size - The angular FWHM assuming a circular cloud, 2A/Area/7r, where the area is 
that of the half maximum value contour shown in the catalog images (Figure [5]). Due 
to our size constraint of 20'along any axis, the selection function for elongated clouds 
drops below 20'. 

• Vlsr - The Local Standard of Rest velocity of the peak of the fit to the average 
spectrum. 

• Vgsr - The corresponding Galactic Standard of Rest velocity defined by 
V G SR=V LS R+220sin(/)cos(6). 

• AV^ - The FWHM in velocity of the fit to the average spectrum. 

• Tpk - The peak brightness temperature after Galactic subtraction. 



2 The online materials are publicly available at 



http://sites.google.com/site/galfahiccc 
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• Nhi - The maximum HI column density calculated using the fit of the spectrum through 
the peak of the cloud in the unconvolved data. Where the fitting failed on the maximum 
spectrum, the average column density is reported as a lower limit and marked with an 
asterisk. 

• S to t - The total flux of the cloud determined by summing over the ROI in the Galactic 
background subtracted data. 

• Population - The classification of the cloud using the scheme detailed in §U 

HVC - (high velocity cloud) |Vlsr| > 90 km s _1 and near a known HVC complex 

GC - (galaxy candidate) |Vlsr| > 90 km s -1 and far from known HVC complex 

CLVC - (cold low-velocity cloud) |V L sr| < 90 km s" 1 and < 15 km s _1 

WLVC - (warm low-velocity cloud) |Vlsr| < 90 km s _1 , AV > 15 km s _1 , and not Q3 

Q3 - (warm, positive low-velocity cloud in the third Galactic quadrant) 
90 >V LS r> km s~\ AV > 15 km s~\ and 270 >/ > 180 



Results 



The conditions for a cloud to be included in this catalog are an angular size smaller 
than 20' and separation from other HI emission both spatially and in velocity. Without 
distance constraints, the clouds identified could ha ye masses that range over m any orders of 
magnitude. From 10" 3 - I P" 1 M m stellar outflows iMatthews et al.l boill l2012h. to 10 1 - 10 2 
M Q recooling 'cloudlets' jHeitsch fc Putmanl |2009j), to 10 2 - 10 4 M clouds at the disc- 
halo interface JFord et al.ll2008l.l201oh. to 10 5 Mp, gas -rich dwarf galaxies (IRyan- Weber et al.l 



20081 ). to > 10 6 M minihalos (IGiovanelli et al.ll2010l ). observations and theories have shown 
that compact neutral hydrogen clouds are a heterogenous group in need of multi-wavelength 
analysis. With this in mind, we have examined the distribution of observable properties and 
separated the catalog into five populations based on velocity, linewidth, and position. We 
recognize that this grouping is inexact and that the populations are continuous across several 
of the boundaries, but we find that there are distinct differences between the populations. 
In this section, we discuss our cloud classifications and describe the observed properties of 
each population. In §[5] we discuss how each population relates to previous observations and 
models. 

Our first step in separating the clouds into populations is to split t he catalog into low- 
and high- velocity clouds at | VlsrI = 90 km s _1 as used in previous studies ( jWakker fc van Woerden 
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19971 ). We chose not to use the deviation velocity (Viev) introduced by IWakkerl (ll99ll ) be- 



cause the GALFA-HI DR1 coverage includes mostly high absolute galactic latitudes and 
the Galactic anti-center where Vlsr— Vd cv - The greatest difference between VLSR,and Vd cv 
in the DR1 region is only 16 km s _1 . The high- velocity clouds are further separated into 
two populations by position and velocity: high-velocity clouds (HVC) that are near known 
HVC complexes as described in §4.11 and galaxy candidates (GC) that are far from known 
HVC complexes. To motivate separating the low-velocity clouds, we note that 83% of the 



neutral medium ( 


Wolfire et al. 


1995 


Putman et al. 


2002: 


Hsu et al. 


2011 



200 31 WNM)and has been observed previo usly (e.g 



de Heii et al.l 120021 : Kalberla fc Haud 



2006|). We use 

this observed boundary to separate the LVCs into cold and warm populations. The cold 
low-velocity clouds (CLVC) have |Vlsr| < 90 km s _1 and AV < 15 km s _1 . The final two 
populations are both warm low-velocity clouds. The warm, low-positive velocity clouds in 
the third Galactic quadrant (90 >V LS r> km s _1 , AV > 15 km s _1 , and 270 > I > 180; Q3) 
are considered a different population from the rest of the warm low-velocity clouds (WLVC) 
due to their distinct velocity distribution (see §4.3. 1[) . 



4.1. High- Velocity Clouds 



We find 692 clouds in our catalog that meet the standard definition of high-velocity 
clouds (HVCs) of | Vlsr | > 90 km s -1 and are not in the galaxy candidates below. HVCs 
are thought to be clouds of gas that reside in the halo of our Galaxy, ranging in dis- 
tance from a few kpc to at least 50 kpc for those associated with the Magellanic System 
(IWakker fc van Woerdenl 119971 : IWakkerl 120041 ) . HVCs are typically found in groupings in 
position-velocity space called complexes. The HVCs detected in this catalog are very small 
compared to the scale of complexes wh ich often span tens of deg rees, and are called ultra- 
compact (size < 20') HVCs fUCHVCs) (brims fc Westmeieril2004f ). 



To determine whether our UCHVCs are associated with known complexes we use the 
updated Wakker & van Woerden (1991; WvW) catalog of HVCs. This catalog has the 
velocities, fluxes, and positions of 616 HVCs across the entire sky. Details on this sample of 
clouds can be found in Wakker (2004). To determine whether the 692 UCHVCs cataloged 
here are associated with known complexes of HVCs, we determine the "distance" in position- 
velocity space between a given catalog UCHVC and a HVC in t he WvW catalog. We 
parameterize distance using the relationship from lPeek et al.l (120081 ): 



D = yje 2 + f 2 {5vf 



(2) 
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where B is the angular distance, 5v is the difference in velocity, and / is a conversion factor 
Peek et al. set equal to 0.5°/km s _1 . By analyzing simulations of HVCs, Peek et al. found 
that clouds from the same complex typically have D < 25°. Figured] is a position-velocity 
diagram where we have colored the HVCs by which complex they are associated with. Most 
notable are the large number of clouds associated with the tip of the Magellanic Stream 
(MS), the anti-center complex (AC), and the WA, WB and WC complexes. Figure [2] shows 
the cloud positions and GALFA-HI DR1 coverage in Galactic coordinates with outlines of 
the HVC complexes. We find that the majority of our HVCs have D < 9° to the nearest 
WvW HVC, with only 35 HVCs having D > 25°. We separate these isolated clouds into a 
sub-population of galaxy candidates discussed in the following section. None of the —400 < 
Vlsr< — 90 km s -1 clouds have D > 35°, leaving 86% of the negative HVC parameter space 
completely empty. 

The HVCs are shown in black in Figures [IJ [7J and |9j Figure [7] shows the linewidth 
distribution for all the clouds in the catalog. We do not separate the HVCs by linewidth 
and note that 83% have AV > 15 km s _1 . The measured properties of the clouds are 
plotted in Figure EJ The size distribution of the HVCs is peaked near 8' with a tail to our 
limit of 20'. The number of HVCs increases at lower peak brightness temperatures, as we 
see for all the populations. There is a sharp decline in clouds below Tpk = 0.3 K, which 
corresponds to ~ 3er in the unconvolved data cubes. The column density distribution peaks 
near ~ 10 19 cm -2 and extends an order of magnitude higher and lower. There is a sharp 
break in the velocity distribution of the HVCs at 90 km s _1 where we define the boundary 
of HVCs and LVCs, indicating some of the LVCs may be low velocity extensions of HVCs. 
We discuss this further in §5.1.11 At the upper ends, the UCHVCs extend to just below -400 
km s _1 (representative of gas at the tail of the Magellanic Stream) and generally under 200 
km s _1 on the positive velocity side. 

4- 1.1. Galaxy Candidates 

There are 27 catalog members which lie at |Vlsr| > 90 km s _1 and are not correlated 
with known galaxies or HVC complexes. These are represented by red in Figures [7J El and 
|9j Of these 27, 8 are at negative velocities and although they lie outside the threshold 
for association discussed in §4.1[ they do appear to be extensions in position-velocity space 
of known HVC complexes. The remaining 19 positive velocity galaxy candidates range in 
velocity from 153 km s _1 to 617 km s _1 , with a median value of 418 km s . A plot of RA 
vs. Vlsr is shown in Figure where this population is shown in red. 

It is likely that many of these high velocity catalog members represent undiscovered 
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galaxies, given their clear bias towards recessional velocities. Supporting this hypothesis 
is their small spatial extent, which has a median size of 5'. While these objects do not 
have peak column densities or brightnesses similar to the 58 known galaxies in our sample 
(average Nhi = 1-69 x 10 20 cm~ 2 and T# = 1.5K), this may be expected for low surface 
brightness dwarf galaxies that have not yet be e n disc overed optically. Many of these are the 
same detections discussed by iGiovanelli et all (120101 ) as possible mini halos. 

We compared the sample of 27 potential galaxies to the complementary Arecibo Legacy 
Fast ALFA Survey (ALFALFA) catalog (Haynes et al. 2011). Much of the DR1 data was 
taken commensally with this survey. Eight of the 27 GALFA-HI galaxy candidates are in 
the ALFALFA coverage region, and 6 of these are in the ALFALFA catalog. For this sample 
of six, Haynes et al. searched for optical counterparts and found 2 to lie outside the clean 
SDSS coverage region, 2 to be of uncertain SDSS association (e.g. displaced in redshift, or 
near objects without known redshifts), and 2 to have no optical counterpart at all in SDSS. 



4.2. Cold Low- Velocity Clouds 



We identify 1245 clouds with |V L sr| < 90 km s" 
temperature and linewidth of a cloud are related by 



T = 21.8K 



AV 

km s~ 



For purely thermal broadening, the 



(3) 



The upper limit on the temperature of a cloud with AV = 15 km s~ x is ~5000K which 
is a typic al lower bound on the temperature of the partially-ionized warm neutral medium 
(WNM) (IWolfire et al.ll2003l ). We separate these low and intermediate velocity clouds into 
750 CNM-dominated cold clouds with average linewidths less than 15 km s~ l and 495 WNM- 
dominated warm clouds with average linewidths greater than 15 km s _1 . As Figure [7] shows, 
the majority of the HVCs i n the catalog are th us classified as warm, as is expected for clouds 
embedded in the hot halo (IWolfire et all 120031 ). 



The cold LVCs are shown in blue in Figures [H |2l [7J and [H From the histograms in 
Figure [9j we see that the cold low- velocity clouds are relatively evenly distributed in velocity 
with 63% having negative velocities. The number of cold LVCs also decreases significantly 
toward ±90 km s _1 , unlike the warm LVCs (as discussed below). Spatially, there appears to 
be an increase in density at a, 5 = 330°, 20°(Z, b ~ 77°, -27°) and an underdensity at positive 
velocities for a = 200 — 250° (general direction of the Galactic center and high positive 
Galactic latitude) (Figure [1]). The cold LVCs are typically smaller than the warm LVCs 
and HVCs with a peak centered at 6'(near our resolution limit), though there is also a tail 
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extending to 20'. This indicates that some of these clouds are unresolved. The linewidths also 
extend to the lower limit of our search space (~3 km s _1 ), indicating there are even colder 
clouds yet to be found. The cold clouds have the shallowest peak brightness distribution 
with the 95 percentile at 1.8K. Similar to the other populations, there is a sharp decline in 
the number of clouds with peak brightnesses below 0.3K, effectively our sensitivity cutoff. 
The steepness of this decline indicates that there may be many more clouds to be found 
with deeper observations. Even though the cold HVCs have higher brightness temperatures, 
they have lower total flux and column density values due to narrower linewidths and smaller 
sizes. 



4.3. Warm Low- Velocity Clouds 

Of the cataloged clouds, 495 have low velocities with |Vlsr| < 90 km s -1 and linewidths 
of AV > 15 km s -1 . As described in §4.2[ this corresponds to temperatures greater than 
5000K. We identify a sub-population of warm LVCs with positive velocities positioned in 
Q3 that we discuss in §4.3.11 and §5.1.11 With the removal of the clouds in Q3, 302 clouds 
remain, 90% of which have Vlsr< km s _1 (magenta in the figures). In Figure [91 we see that 
the velocity distribution of the warm LVCs continues to -90 km s" 1 where we separate them 
from the HVCs, indicating that these two populations are most likely blended in velocity 
or even partially part of the same population. There is a sharp drop in clouds approaching 
km s" 1 due to confusion with larger structures of low velocity gas. The angular size 
distribution is similar to the HVCs with a median value of 9' and extending to the limits 
of the catalog. The total flux and column density distributions are also similar to the HVC 
distributions. We observe a steeper peak brightness distribution than for the cold LVCs with 
95% of warm LVCs below 1.0K. The velocity widths show a sharp cut-off at our break at 15 
km s _1 indicating there is most likely a continuum of clouds in temperature. There appears 
to be overdensities of warm clouds located at 200°-250° and 330° in RA, the locations of 



IVC complexes K and the PP Arch, respectively (jWakkerll200ll ). 



4-3.1. Warm, Positive- Velocity Q3 Clouds 

The vast majority of warm clouds with 90 >Vlsr> km s -1 are located in the third 
Galactic quadrant (RA=5-14hr; green in the figures). Since our coverage does not include 
the fourth quadrant (see Figure [2] for the GALFA-HI DR1 coverage in Galactic coordinates), 
the third quadrant is the only area of the sky where infalling (or static) clouds may have 
positive LSR velocities due to Galactic rotation. These clouds are below the traditional 
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velocity cut for HVCs, but otherwise have similar properties to the HVCs. In Figure [U we 
plot the velocities that correspond to Vqsr = km s _1 to illustrate the effect of Galactic 
rotation. In the Galactic Standard of Rest (GSR), these clouds have velocities near or below 
zero, as do the HVCs associated with the WA, WB, and WC Complexes. This suggests that 
the positive velocities in the LSR frame are a product of Galactic rotation, and that these 
clouds have low or negative velocities relative to the Galacti c plane. In po sition, the warm 



Q3 clouds appear to be an extension of the IV-WA complex (lWakkerll200ll ) which is located 



at similar Galactic latitudes but in the fourth Galactic quadrant. 

From Figure |9l the warm Q3 clouds have size, total flux, and column density distri- 
butions similar to the rest of the warm LVCs and the HVCs. The linewidth distribution is 
similar to the warm clouds with a sharp cut at 15 km s _1 and extending beyond 30 km s _1 . 
The warm Q3 clouds are the weakest population in brightness temperature with a 95 per- 
centile of 0.8 K. The velocity distribution shows that the decrease in the number of clouds 
toward Vlsr=0 km s _1 begins at a higher velocity than the rest of the warm or cold cloud 
populations. This may be partially due to only including a limited range of positions for 
this population in the direction of positive Galactic rotation; however, there does appear to 
be a larger gap between the warm Q3 clouds and the Galactic disc in Figure [T] than the rest 
of the warm LVC population. 



5. Discussion 

In this section we discuss how the different populations of clouds we have defined in 
fit into current theories and how they relate to previous observations. 



5.1. UCHVCs 

UCHVCs are especially of interest because it has been suggested that small HVCs 
may represent phy sically large structures at megaparsec distanc es, much farther than the 



classical complexes (IBraun fc Burton! Il999t iGiovanelli et al.ll2010l ). It is very hard to obtain 



direct distance information for such small clouds, given the unlikely chance of an overlap 
with a halo star; however, if the UCHVCs reside near larger complexes in position-velocity 
space, we expect they are associated with these relatively nearby structures, rather than 
being independent clouds at much larger distances. 

As discussed in §E1 all but a few UCHVCs have D < 25° to an HVC in the WvW 
catalog. At negative velocities, this leaves 86% of the available position-velocity space with 
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-400 <V LS r< -90 km s" 1 empty. The presence of UCHVCs found only near other known 



Binnev et al. 


2009; 


Joung et al. 


2011) 



thermally unstable to linear perturbations, we would expect small neutral condensations to 
appear all over phase space, rather than only concentrated toward known, larger HVCs. This 
result instead supports the idea that some HVCs may be seeded by larger perturbations in 
the halo, such as filamentary streams of gas impinging on the Galaxy fro m the intergalactic 
medium or satellite galaxies ( jKeres fc Hernquistll2009t IJoung et al.ll201ll ). 



If UCHVCs are indeed associated with larger clouds, they are most likely at similar 
distances. The median UCHVC properties indicate a mass of ~ 2OOM and a physical size 
of ~ lOpc at a distance of lOkpc and a mass of ~ 2 x 10 4 M Q and a physical size of ~ lOOpc 
at lOOkpc. These distances were chosen based on the distance constraints to many HVC 
complexes ( iThom et al.l l2006t IWakkerl l200ll ; IWakker et al. l2008h . and the likely distance o f 



the tail of the Magellanic Stream given simulations (jBesla et al.ll2010t IConnors et al.ll2006l ). 
The mass of these clouds indicate they are not likely to survive t he trip to the Galactic dis k 
and will instead become part of the multi-phase Galactic halo ( jHeitsch fc Putmanl 120091 ). 
The exception would be if some of the UCHVCs are embedded within extended, highly 
ionized shells. 



5.1.1. Warm, Positive-Velocity Q3 Clouds 

In §4.3.11 we state that the vast majority of the positive velocity warm low velocity 
clouds (0 <Vlsr< 90 km s _1 ) are located in the third Galactic quadrant and are likely to 
be associated with the W HVC complexes. In Figure [T] we have plotted the Vlsr velocities 
corresponding to Vqsr= km s _1 as a curved line. Clouds on this line could be at rest 
relative to Galactic rotation. The majority of the positive-velocity warm clouds in Q3 are 
bounded by the Vgsr= km s -1 line indicating these clouds may be infalling clouds with 
positive LSR velocities due to Galactic rotation. 

The HVC complexes WA, WB, and WC (e.g. WvW) are in the same region of sky 
as the warm Q3 clouds, and if we extend the associations with HVCs from ^4.11 below 90 
km s" 1 , many of the Q3 clouds would be considered part of these complexes. iThom et al. 



( 120061 ) measured a direct distance constraint to clouds in the W complexes and found a 
distance of ~ 9 kpc. This results in a z-height of 7 kpc and galactocentric distance of 12 kpc. 
We consider three scenarios for a cloud at this distance at (a,5) = (150°,15°): a non-rotating 
cloud, a co-rotating cloud, and a co-rotating cloud with a vertical lag. For a non-rotating 
cloud (Vgsr=0), we would observe Vlsr— 95 km s _1 due to the rotation of the Galaxy. If the 
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cloud were co-rotating at 220 km s _1 , we would observe Vlsr— 10 km s _1 . If the cloud were 
co- rotating, but with a vertical lag of 20 km s -1 kpc -1 , we would observe Vlsr— 65 km s" 1 . 
All of these calculations are without any vertical infall which would lower the observed Vlsr- 
The velocities we observe for the warm Q3 clouds are consistent with non-rotating clouds 
with infalling velocities of ~10-70 km s" 1 or co-rotating clouds in a laggi ng halo. The warm 
Q3 clouds may fall under the category of the 'low- velocity halo clouds' (IPeek et al.ll2009l ). 
and may also represent a bridge connecting the HVC W complexes to the Galactic disc. 



5.2. Cold Low- Velocity Clouds 

The cold, low-velocity clouds are certainly associated with the disc of the Galaxy. There 
are greater numbers of these clouds at lower velocities and there is no apparent relationship 
to the GSR frame lik e for the warm Q3 cloud s ( §5.1.ip . The cold LVCs may be related to 



the clouds studied by iFord et all (120081 . 120101 ) who identified a population of discrete, cold 



clouds 600-1000pc out of the plane of the Galaxy that are co-rotating with the disc. Their 
measurements were made on a set of clouds with Galactic latitude less than 20° and at a 
distance of ~8kpc in the first and third Galactic quadrants. Ford et al. determined distances 
by selecting clouds above the tangent point velocity, which is the maximum velocity allowed 
by Galactic rotation toward sightlines of the inner Galactic plane. The GALFA-HI DR1 does 
not include the inner Galactic plane so we cannot use the tangent point method to determine 
distances. Much of this catalog is at high Galactic latitudes where the Ford clouds would be 
closer (< lkpc) and therefore appear substantially larger, up to two degrees, though many 
of their clouds were unresolved. The masses of the Ford sample range from ~ 100 to ~ 5000 
M Q , while the cold LVCs from our catalog would have masses of ~ 1 M Q at a similar vertical 
position. High latitude analogs of the Ford clouds would not be included in this catalog due 
our size limit of 20', while analogs to the clouds in this catalog would not have been detected 
by Ford et al. due to the small sizes and masses. However, the clouds in the Ford sample 
have a median linewidth of 10.6 km s -1 , while our cold clouds have a median linewidth of 8.4 
km s _1 , indicating both populations have typical temperatures less than 5000K. The cold 
LVCs in this catalog are not physically the same as the Ford clouds, but may have similar 
temperature, vertical structure, and kinematics. 

Two possible origins of the cold LVCs are a Galactic fountain process or recooling 
disc-halo interface gas. In a Galactic fountain model, gas is expelled out of the plane of the 



Gala xy by some energetic process and then rains back down out of the halo (jShapiro &: Field 



19761 ) . Whether the ejected gas is localized to regions of star formation, or rapidly blends 



with existing warm/hot lower halo gas and condenses later is still unknown. This will be 
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investigated further with models and a larger popu lation of GALFA-HI clouds in the future. 
Regarding recooling gas, iHeitsch fc Putmanl (120091 ) observed small clouds forming from the 
warm remnants of disrupted HVCs in their simulations as the clouds became buoyant close 
to the disk. The large number of positive-velocity cold LVCs requires some process to drive 
them away from the disc, be it the buoyancy of accreting gas or feedback from star formation 
in the disc. 



5.3. Warm Low- Velocity Clouds 

After classifying the positive-velocity warm clouds in the third Galactic quadrant as 
possible low- velocity halo clouds, we see that 90% of the remaining warm low- velocity clouds 
have negative velocities. This agrees well with previous studies of intermediate velocity 
(|Vlsr| = 40 — 90 km s~ x ) clouds that found most o f this gas is at negative velocities and has 



a vertical distance ranging from 0.5 to 3 kpc (See I Albert fc Danlyl (120041 ) for an extensive 
discussion). The warm low- velocity clouds detected here may be associated with IVCs. If 
these clouds are within 3kpc of the disc, they would have masses of ~ 1 — 100M Q . The 
origin of the clouds is uncertain but the higher metallicities for some of the large IVCs 



(IWakkerll200ll ) suggests a relation to processes within the Galactic disc. As shown in Figure 
O the velocity and line width distributions are continuous across the 90 km s -1 and 15 
km s _1 boundaries, respectively. This suggests that the warm LVCs may be partially low- 
velocity HVCs and partially warm or turbulent cold LVCs. This is different from the cold 
LVCs that appear separate from the HVCs in velocity. From this reasoning, some of the 
warm LVCs are most likely the remnants of disrupted HVCs, while other warm LVCs may 
have come from, or have yet to become, cold LVCs. In any case, the negative velocities of 
these clouds indicate they will be future star formation fuel. 



6. Summary 

Using a custom cloud detection algorithm, we have identified 1964 compact (< 20') neu- 
tral hydrogen clouds with velocities typically between |Vlsr| = 20 — 450 km s _1 , linewidths 
of 2.5 — 35 km s _1 , and column densities ranging from 1 — 35 x 10 18 cm -2 . At high velocities, 
we observe the vast majority of clouds are associated with previously known HVC complexes. 
At low velocities, we observe populations of cold and warm clouds that appear associated 
with the disc of the Galaxy, as well as a group of warm clouds at positive velocities that may 
be low-velocity halo clouds associated with HVC complexes WA, WB, and WC. The catalog 
extends down to our limits spatially, spectrally, and in sensitivity, indicating that while we 
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have discovered a large number of new small structures, there are still smaller, colder and 
fainter clouds to be detected. 

Although our search criterion was simply for clouds smaller than 20' in the GALFA-HI 
DR1 area, we detected clouds in vastly different environments. From potentially undiscovered 
galaxies, to HVCs associated with the Magellanic Stream at 100 kpc, to complexes of clouds 
at < lOkpc, to clouds associated with the Galactic disc, we find differences in the distribution 
of clouds indicating a range of origins. While the HVCs associated with complexes may be 
sheared off the larger clouds by instabilities, the LVCs could have been ejected from the disc 
and some may represent a link to HVCs. The distances to the clouds in this catalog may 
cover several orders of magnitude, so the masses may range from less than a Solar mass for 
clouds within the Galactic disc, to greater than 10 4 M for HVCs at the tip of the Magellanic 
Stream. 

This catalog is ripe for mult i- wavelength analysis and follow up observations. Finding 
optical counterparts to the galaxy candidates will provide distances and masses, and work is 
already underway to identify potential new ultrafaint dwarf galaxies with the catalog (Grce- 
vich et al., in prep.). Infrared data will allow us to determine the dust content of these clouds 
and may disentangle the halo clouds from the disc clouds. We are also pursuing Herschel 
spectral line observations to investigate the possibility of cooling within the clouds, and will 
look for any fortunate overlap with stars to provide direct distance and potentially metallic- 
ity constraints. Finally, future GALFA-HI observations combined with Galaxy models will 
further constrain the origin and role of the smallest and coldest population of individual HI 
clouds yet found. 
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0707679, 0709347, 0917810, NSF CAREER grant 094059, NASA Herschel grant 1369759, 
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Table 1. DR1 GALFA-HI Cloud Catalog 



Cloud 


RA 


DEC 


1 


b 


"Max 


Size 


Vlsr 


Vgsr 


AV 


rpPk 

l b 




Stot 


Population 


ID 


(h:m:s) 


(d:m) 


(deg) 


(deg) 




(arcmin) 


(km s' 1 ) 


(km s _1 ) 


(km s- 1 ) 


(K) 


(10 18 cm- 2 ) 


(Jy km s- 1 ) 




001.4+35.7+030 


00:05:29 


35:40 


112.54 


-26.27 


16.4 


6.5 


30.98 


213.19 


8.42 


1.29 


16 


12.2 


CLVC 


001.9+16.3-350 


00:07:37 


16:17 


107.91 


-45.30 


14.8 


8.1 


-346.88 


-199.62 


37.42 


0.45 


32 


13.8 


HVC 


096.3+35.2-079 


06:25:13 


35:09 


178.59 


10.34 


22.7 


12.2 


-80.00 


-74.66 


4.26 


1.69 


10 


8.5 


CLVC 


109.7+14.7+063 


07:18:49 


14:39 


202.74 


12.60 


50.7 


19.5 


64.70 


-18.31 


21.08 


0.83 


33 


15.2 


Q3 


143.1+27.1-119 


09:32:29 


27:06 


201.03 


46.04 


13.9 


7.0 


-119.42 


-174.21 


4.78 


0.48 


3 


2.7 


HVC 


175.8+32.2-070 


11:43:21 


32:14 


189.75 


74.38 


8.3 


4.7 


-70.12 


-80.15 


3.10 


0.34 


2 


2.0 


CLVC 



Note. — An excerpt from the full catalog table correspond to the clouds in Figure[5] The velocity in the cloud ID is determined from the peak position in the 
convolved cube and varies slightly from Vlsr which is the velocity from the Gaussian fit. The full table is available at jhttp: //s ites. google.com/site/galfa hicccj 
*Average column density is reported where the fit through the maximum of the cloud fails to converge. 
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Table 2. Median Properties of Cloud Populations 



Parameter 


HVC 


Gal 


Cold 


Warm 


Q3 


V L sR.(km s" 1 ) 


-166 


266 


-38 


-59 


61 


T pk (K) 


0.4 


0.4 


0.5 


0.4 


0.4 


Av (km s _1 ) 


22.0 


21.3 


8.4 


20.0 


20.0 


N H i (xlO 18 cm~ 2 ) 


11 


8 


6 


9 


9 


Angular Size (') 


9 


5 


7 


9 


9 



Table 3. 90% Range of Properties of Cloud Populations 



Parameter 


HVC 


Gal 


Cold 


Warm 


Q3 


V L sR.(km s" 1 ) 


-389 ->■ 126 


-237 ^ 612 


-79 ->■ 64 


-89 ->■ -48 


21 ->■ 86 


T pk (K) 


0.2 ->■ 1.1 


0.3 ->■ 1.1 


0.3 ->■ 1.8 


0.2 ->■ 1.0 


0.2 ->■ 0.8 


Av (km s _1 ) 


6.3 ->■ 34.3 


2.5 ->■ 69.8 


2.5 ->■ 14.2 


15.4 ->• 33.8 


15.4 ->■ 32.8 


N HI (xlO 18 cm" 2 ) 


2 ^ 34 


2 ^ 70 


1 ->• 19 


2 ^ 28 


3 ^ 26 


Angular Size (') 


5 ^ 15 


3^8 


4 ^ 12 


5 ^ 14 


5 ^ 14 
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Fig. 1. — Position-position- velocity distribution of clouds with Vlsr< ±250km s _1 . The 
cloud populations are colored as in Figure [7J Bottom Panel: Position-position plot showing 
the GALFA-HI DR1 coverage. The hashed region is not included in DR1. The grey hor- 
izontal lines indicate the breaks between the declination regions used in the upper panels. 
Upper Panels: Position-velocity plots showing clouds in the declination range indicated. The 
shading indicate s the relative amount of observations in the region. The greyscale image is 
the LAB survey (IKalberla et al.ll2005l ) integrated along the declination range indicated. The 
curved line follows the velocities corresponding to Vgsr= km s -1 , where a cloud at rest 
compared to Galactic rotation would lie. 
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Fig. 2. — A Hammer-Aitoff plot of the positions of all of the cataloged clouds in Galactic 
coordinates. Cloud type is indicated by color, with black for HVCs, red for galaxy candidates, 
blue for cold LVCs, pink for warm LVCs (excluding the Q3 clouds), and green for warm, 
positive-velocity Q3 clouds. The coverage of the survey is indicated in grayscale, with the 
darkest gray indicating > 50 seconds per beam of integration time. High velocity cloud 
complexes are shown in outline, with the W complexes shown in green, complex C shown in 
blue, and the Magellanic Stream shown in orange. Other HVCs are shown in black. 
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Fig. 3. — An example of the Difference of Gaussians kernels used in the convolution step of 
our algorithm. This is a two dimensional slice through one of the three dimensional kernels. 
Here, the positive Gaussian has a FWHM of 7', while the negative Gaussian has a FWHM 
of 8.4'. The Gaussians are scaled such that the integral of the 3D kernel is zero. 
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Fig. 4. — The sensitivity to recovering clouds of given sizes and velocity widths given our cho- 
sen optimal kernels. The cut-off at the lower end in size and velocity is set by the resolution 
of the data. Excluding the clouds that are both largest (20') and widest (20 km s -1 ), where 
the sensitivity drops by 70%, the sensitivity drops to no less than 80% in the parameter 
space shown. 
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10 20 30 40 50 -400 -380 -360 -340 -320 -300 -90 -85 -80 -75 -70 

Velocity (km/s) Velocity (km/s) Velocity (km/s) 
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40 60 80 100 -130 -125 -120 -115 -110 -105 -75 -70 -65 

Velocity (km/s) Velocity (km/s) Velocity (km/s) 



Fig. 5. — A random sample of the DR1 GALFA-HI cloud images and spectra (the full version 
is available online). The image is produced by scaling each velocity channel by the result 
of the fit to the average spectrum (upper spectrum) and then integrating along velocity. 
This de-emphasizes nearby emission at velocities away from the peak of the cloud while still 
increasing the fidelity of the image over an unintegrated slice. The contour follows the half- 
maximum of contiguous emission surrounding the peak position and is the reported area of 
the cloud. The upper spectrum is the average over the Region of Interest with the four- 
parameter Gaussian fit (see §3.11 and §3.4p . The lower spectrum is the spectrum through 
the peak in the unconvolved data. If a Gaussian fit was successful on this spectrum, it is 
plotted. The full catalog is available at |http: / /sites.google.com/site/galfahiccc| 
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Fig. 6. — The ratio of our sensitivity function, S (T pk , a, b/a, 5v, t), to the signal detected in 
the convolved data, a c i, as a function of various parameters of the injected clouds. We find 
we can accurately predict a c i using Equation CD 
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Fig. 7. — Linewidth and VLSRdistributions for all clouds. The differently colored symbols 
correspond to the five populations - HVCs (black triangles), galaxy candidates (red dia- 
monds), cold LVCs (blue Xs), warm LVCs (pink open squares), and warm Q3 LVCs (green 
filled squares). This plot best illustrates where the populations are separated. See Figured] 
for the velocity and linewidth distribution for each population. 
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250 200 150 100 50 

RA [degrees] 

Fig. 8. — We show the distribution of all cataloged clouds in RA vs. Vlsr, highlighting the 
HVCs. Clouds the do not meet the simple HVC criterion of |Vlsr| > 90 km s _1 are marked 
in gray. Clouds with a D (see Equation [2]) to the closest WvW catalog HVC less than 25°are 
marked in the color denoting that nearest complex or grouping in the color bar at right. 
Clouds marked in red have D greater than 25° to any HVC in the WvW catalog. Note the 
strong asymmetry across Vlsr= 0km s _1 in unassociated clouds. 
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Fig. 9. — The distributions of angular size, peak brightness temperatures (T p k), column den- 
sity, LSR velocity, total flux and linewidth for all cataloged clouds, separated by population 
- HVCs (black, lower pane), galaxy candidate (red, second pane), cold LVCs (blue, center 
pane), warm LVCs (pink, fourth pane), and warm Q3 LVCs (green, upper pane). To aid in 
comparing the different populations, each histogram is scaled so that the fractional number 
of clouds per bin is consistent for each plot. The y-axis values are the mid- values for each 
histogram in number of clouds. 



